Steady-State Chemotactic Response in E. coli 



O 

o 

> 
O 

ON 

(N 



u: 
6 : 
X) : 

(N : 
> ■ 

: 

(N . 

\o • 

d ■ 

: 

o ■ 



X 



Yariv Kafri 1 and Rava Azeredo da Silveira 2 
1 Department of Physics, Technion, Haifa 32000, Israel 
2 Department of Physics and Department of Cognitive Studies, 
Ecole Normale Superieure, 24, rue Lhomond, 15005 Paris, France 

The bacterium E. coli maneuvers itself to regions with high chemoattractant concentrations by 
performing two stereotypical moves: 'runs', in which it moves in near straight lines, and 'tumbles', in 
which it does not advance but changes direction randomly. The duration of each move is stochastic 
and depends upon the chemoattractant concentration experienced in the recent past. We relate 
this stochastic behavior to the steady-state density of a bacterium population, and we derive the 
latter as a function of chemoattractant concentration. In contrast to earlier treatments, here we 
account for the effects of temporal correlations and variable tumbling durations. A range of behaviors 
obtains, that depends subtly upon several aspects of the system — memory, correlation, and tumbling 
stochasticity in particular. 



Chemotaxis refers to directed motion in response to 
chemical signals and has been extensively studied in the 
bacterium Escherichia coli (E. coli) E. coli is con- 
fined to two stereotypical moves. When its flagellum mo- 
tors turn counterclockwise (looking at the bacteria from 
the back), the bacterium moves in near-straight lines 
termed 'runs' whose direction is limited by rotational dif- 
fusion. This motion is interrupted by periods of 'tumble' 
which occur when the motors turn clockwise: the bac- 
terium does not translate but instead rotates about itself 
in a random fashion, and thus reinitializes the direction 
of the next run. Some amount of correlation between 
successive run directions yields an average angle shift of 
68° [H , as compared with the 90° in an uncorrelated case. 
Tumble durations are short, of the order of 0.1 s, with 
respect to runs which last about 1 s @, 0] • 

Bacteria modulate their whereabouts in response to 
their chemical environment. The small size (~ 2/um) of 
E. coli rules out sensing spatial gradients: in the time it 
takes the bacterium to move by its own size, chemicals 
diffuse in a region ten times larger Q • Instead, E. coli 
relies upon temporal integration: it calculates a spatial 
gradient by integrating the concentration of chemicals 
over its recent history; it then uses the resulting quantity 
to modulate run and tumble durations. Much work has 
focused on this 'algorithm', namely, on the filter of tem- 
poral integration and on its relation to the probability of 
run or tumble. P, 0, Q . 

With knowledge of this stochastic algorithm, one 
would like to predict the distribution of trajectories of 
a bacterium or, equivalently, the behavior of a (non- 
interacting) population. And, in particular, one would 
like to elucidate which aspects of the single-bacterium al- 
gorithm ensure population performance. Here, we focus 
on the steady state and we ask the following questions. 
Given a chemoattractant (or chemorepellant) concentra- 
tion and a single-bacterium stochastic algorithm, what is 
the shape of the steady-state population density? How 
does it depend upon the details of the single-bacterium 
system and which of these are qualitatively relevant? 



Because of the single-bacterium stochasticity, the prob- 
lem may be viewed as a biased random walk problem. 
The memory involved in temporal integration and the 
variable tumble duration, however, make the problem 
more difficult and more interesting. In particular, these 
induce correlations between run duration and local bac- 
terium density at the run's starting point, which we took 
into account. All these effects yield a rich macroscopic 
behavior in the steady state that depends subtly upon 
the form of the single-bacterium response filter. In par- 
ticular, (i) the usual bi-lobe filters that turn temporal in- 
tegration into spatial comparison may or may not lead to 
accumulation in favorable regions, depending upon their 
shape and the interplay of time scales 0, (ii) cor- 
relations result in a non-local dependence of the prob- 
ability density upon the environment, due to the pres- 
ence of memory in the dynamics; (Hi) when tumble is 
non-instantaneous, bacteria may aggregate in favorable 
regions in their tumbling phase. Suprisingly, this last ef- 
fect occurs even for filters that are purely local in time. 
Our results are derived in one spatial dimension, as in 
Refs. 7, 8], and fodder a long-standing debate [J la, H, H| ; 
specifically, correlations and tumble duration variability 
had been neglected in earlier treatments 0, H, H, 0, @, El • 

E. coli climbs up chemical gradients by modulating 
run and tumble durations as a function of chemoattrac- 
tant concentration, c dUl. (Henceforth, we use the term 
'chemoattractant' indifferently to refer to both chemoat- 
tractant and chemorepellcnt. Below, we discuss the dif- 
ferences in responses to 'positive' and 'negative' stimuli.) 
Run durations are Poissonian, with probability 
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to switch from run to tumble between times t and t + dt 
Q . Here, T [c] is a functional of the chemical concentra- 
tion, c(t'), experienced by the bacterium at times t' < t; 
it results from a linear temporal filtering followed by a 
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static rectification non-linearity, as 

F[c]=cf>(^J dt'R(t-t')c(t')j , (2) 

where the functions <fi (•) and R(t) summarize the action 
of the biochemical machinery that processes input signals 
from the environment []]]. If <j> (•) is non-singular, it may 
be linearized, as 

^Un[c] = [ dt'R(t-t')c{t'), (3) 

J — oo 

where an additive constant is absorbed in a redefinition of 
To (in Eq. (fT])) and a multiplicative constant is absorbed 
in a rescaling of R(t). Experimental work [3J] suggests 
instead a thresholding non-linearity [9( , well fitted by the 
form 

•fnlin [C] = [Hn [c]] + , (4) 

where 

[x] + ={!!x>r (5) 

The response filter, R(t — t'), was measured in clas- 
sic experiments on wild-type bacteria, in which puffs of 
chemoattractant were presented to a single bacterium, 
effectively replacing c (t 1 ) by a delta- function which al- 
lowed one to resolve for R (t - t') 0, i]. These exper- 
iments yielded a bimodal shape for R(t — t'), with a 
positive peak around t' ~ 0.5 sec and a negative peak 
around t' ~ 1.5 sec. The negative lobe is shallower than 
the positive one and extends up to t' ~ 4 sec, beyond 
which it vanishes and to a good approximation satisfies 
J R(t')dt' — 0. The estimated value of To is about 1 sec 
(see Fig. 1 for an illustration). 

Tumble duration also is modulated stochastically, in 
close analogy to run duration behavior Earlier theo- 
retical work has mostly treated tumble as instantaneous 
0, H, [9j] . We treat tumble duration as a Poisson variable 
with rate 1/tt but, for the sake of simplicity, we ignore 
any dependence of the latter upon the chemical environ- 
ment. While this comes short of a full description, the 
mere allocation of a non-vanishing duration to tumble 
brings in qualitative consequences, as discussed below. 

The bi-lobe shape of the response filter points to a 
simple mechanism: it enables the bacterium to perform 
a coarse-grained temporal derivative of the chemical con- 
centration it experiences. If the gradient is positive, then 
the run duration tends to increase; if the gradient is 
negative, then the run duration tends to decrease (in 
the linear case of Eq. ([3])) or is unmodulated (in the 
threshold-linear case of Eq. (j4|)). However, the connec- 
tion between simple arguments such as this and quan- 
titative results is far from immediate. Reference |7| ar- 
gues that a single-lobe, even punctual temporal filter, 



as R(t — t') = x$(t — t') with \ positive, leads to a 
net bias toward increasing chemoattractant concentra- 
tion. In fact, the analysis suggests that the response is 
strongest if the filter is local in time, with t' = 0, and 
that a delayed response (t' > 0) or any addition of a neg- 
ative contribution, akin to the bi-lobe shape measured 
experimentally, weakens the bias. The arguments devel- 
oped in Ref. [7j concern the instantaneous dynamics of 
a bacterium and ignore the spatially varying buildup of 
probability in time; they apply, for example, to a tran- 
sient situation in which the probability density is flat. 
Reference 0| contrasts transient and steady-state behav- 
ior and argues that while a positive filter is most favor- 
able for climbing chemical gradients in an initial transient 
phase, a negative filter is favorable for steady-state accu- 
mulation in advantageous regions. Finally, it argues that 
the typical bi-lobe shape of the linear filter, R(t), may 
derive from a constrained optimization involving tran- 
sience and steady state. While Refs. 0] present a 
number of interesting ideas and go some length into ex- 
plaining chemotaxis statistically, they make a number of 
limiting assumptions. First, they disregard the correla- 
tion between run duration and probability density in a 
given region. Second, they assume instantaneous tum- 
ble. Third, the constrained optimization in Ref. [H is 
somewhat ad hoc. 

In the remainder of the paper, we proceed as follows. 
First, we write equations that govern the steady-state 
density of (non-interacting) bacteria (equivalently, the 
bacterium probability density); second, we derive the lat- 
ter analytically in the linear model (Eq. ||3J)) and numer- 
ically in the non-linear model (Eq. ((4])). We are after the 
density 

N(x)=N R (x) + N T (x) , (6) 

where Ni(x)dx is the number of bacteria lying between 
x and x + dx in the steady-state, and the subscripts R 
and T refer to run and tumble respectively. As a natu- 
ral way to incorporate correlation, we borrow four inter- 
mediate quantities: n r ~^ R {x)dx, the number of bacteria 
that switch from tumble to rightward run between x and 
x + dx per unit time; n' F ^ R (x)dx, the number of bacteria 
that switch from tumble to leftward run between x and 
x + dx per unit time; n R ^ T (x)dx, the number of bacteria 
that switch from rightward run to tumble between x and 
x + dx per unit time; n R ~* T (x)dx, the number of bacteria 
that switch from leftward run to tumble between x and 
x + dx per unit time. The rightward and leftward fluxes 
are given by 

d xJ+ (x) = n T + - R (x)-n^ T (x), (7) 
d x j-(x) = n T - B (x)-n R - T (x) . (8) 

If N + (x) and N-(x) are the densities of rightward and 
leftward running bacteria respectively, then j+(x) = 
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vN + (x) and j-(x) — —vN_(x), and the total density 
of running bacteria, Nr(x), obeys 

d x N R (x) = d x [N+(x)+N-(x)} (9) 
= -[nl^ R {x)-n^ T {x) 



,T^R 



(x)+n R ^ T (x)] 



(10) 



Tumbling bacteria retain some memory of their recent 
run direction; we call q the probability that a tumble 
causes a run direction change, and treat it as a parame- 



ter in our model. Thus, n^T >R {x) = (1 — q) ; 



,R-*T 



(X) 



qn R r- T (x) andn 1 ^ R (x) = qn R ~^ T (x) + (1 — q) n R ~ >T {x), 
so that Eq. dTUJ) simplifies into 

d x N R (x) = * [n R - T (x) - n R - T (x)] . (11) 



Within our assumption of unmodulated tumble rate, in 
the steady state the density of tumbling bacteria, Nt(x), 
reads 



Nt(x) — tt [r 



R^T 



(x) 



,R-^T 



(*)] 



(12) 



As a final simplifying assumption, we posit that mem- 
ory is erased at tumble-to-run switches. This assumption 
may not be validated by data -3] , but it is unclear whether 
it improves or suppresses chemotaxis with respect the no- 
erasure case. Equation ([3]) becomes 



^iin [c] = / dt'R(t 
J t 



t'W), 



(13) 



where to is the time of last switch, and the run-to-tumblc 
switch probability, dt/T(t,to), is now a function of both 
t and to- Alternatively, this probability can be expressed 
in terms of the initial and final positions of the run, y 
and x respectively, as dx/vT(x,y). We now have all the 
elements in hand to write the steady-state equations that 
govern density and keep track of correlations, as 



x) 



n R ^ T {x) 



dyn 1 + - R (y)p + (x,y), (14) 

Xj 

dynZ-* R (y)p-{x,y); (15) 



these express the fact that tumbling bacteria result from 
running bacteria that switch to tumbling mode. Here, 
p+(x,y)dx and p^(x,y)dx are probabilities that a bac- 
terium, which tumbled last at y, tumbles again between 
x and x + dx (and not before), for x > y and x < y 
respectively. These probabilities are given by 



P±(x,y)dx = exp dy' — ^ -^ j 



dx 



(16) 



We choose to illustrate our results with steps of 
chemoattractant concentration, c(x) = £8(x) (£ > 0), 



where 9(x) denotes the Heaviside function. In the lin- 
ear model, it is handy to focus upon singular response 
functions with R&{t) cx 5(t — A/v) , or equivalently 
in space coordinates, Ra(%) — Xa°~{x T A) (with a 
minus (plus) sign for rightward (leftward) runs). One 
can then derive solutions for more general cases as lin- 
ear superpositions of solution for singular response func- 
tions. We treat the linear model perturbatively in the 
strength of bacterium response; specifically, we assume 
a regime with qa = £xa *C 1. (For otA — 0, there 
is no chemotaxis.) Expanding Eq. (| 16[) to first order 
in a a, we solve the steady-state Eqs. pT|) and (fl5|) 
for the intermediate quantities n R ^ T . From these, we 
derive the incremental running and tumbling bacterium 
densities compared to the densities far to the left of the 
chemoattractant step: 5N^(x) = N^(x) — N^(~oo) and 
5Nft(x) = Nft{x) - Nfi(-ac). 

Because of the singular response function and the dis- 
continuity in chemoattractant density at x = , our so- 
lutions have singular points at x — ±A. We find, for 
x < -A, 



5N£(x) = SN^(x) = 0; 



for -A < x < A, 

5N${x) = -2aq 



a A {x + A) _ A 



/VT 



V 2 T 



(17) 



(18) 



A, ., _ .."AT / | ,,„(A + x)\ ^-A/vro.^Qj 



5N£(x) = -a^^ [l + 2q 



VT 



VT 



for x > A, 



5N*(x) 



6Nfi(x) 



~4aq 



~ — ' 



-A/fro 



-a 



oiatt 



1 



, A 
4q | e 

ur , 



A/vtq 



VTo 

+ 2^Aj ^ SNr W 



(20) 



(21) 



here a is positive a constant that sets the overall density 
of bacteria. From Eq. (|2H)l . running bacteria accumulate 
to the right if «a < , as long as the 'response memory' 
is non- vanishing (A ^ 0). Accumulation is strongest for 
A = vto, i- e - i when the response memory, A/v, is 
comparable to the typical run duration, r . We note also 
that accumulation vanishes if q = 0; indeed, in this case 
bacteria do not change their run direction after tumble 
and, hence, behave roughly as if there were no tumbles 
whatsoever. As typically tt <C to (experimentally, for 
E. coli, tt ~ to/10), Eq. (|2~Tj) implies that SNft is dom- 
inated by SNft. However, the reverse occurs in the par- 
ticular case with small response memory A/v < tt/2, 
i. e., when the typical tumble duration exceeds the re- 
sponse memory. In this case, bacteria may accumulate 
to the right (if «a < 0) even for a response function 
purely local in time (with A = 0) — a possibility over- 
looked in earlier work that treat tumble as instantaneous. 
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In this tumbling-dominated regime, bacteria accumulate 
at favorable tumbling sites while the uniformly populated 
runs serve as a way to explore potentially favorable tum- 
bling positions. 
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FIG. 1: Numerical results in the non- linear model with 
c(x) = W~ s 8(x). All quantities are given in arbitrary units. 
To obtain 6Nr(x) and 5N T {x), Eqs. (|mi5JI . for the non- 
Unear model, were solved iteratively on a computer, using 
a discrete lattice of size 800, with v — Wfim, q = 0.4, 
tt = 10/17s and dx = 0.05/im. The results were uni- 
formly rescaled for convenience and we display the region 
in which the bacterium density varies, about x = 0, the lo- 
cation of the chemoattractant step. Three different values 
of ro (in seconds) are indicated in the figure. The corre- 
spondence between the bacterium density and the value of 
ro is indicated by the solid and open symbols, in both fig- 
ure and inset in which the response function is illustrated. 
The functional form of the response function was chosen as 
R(t) = (240texp(-200t/17) - 29 At exp(-7£/17))/17 which 
satisfies J Q R(t)dt — 0. 

Our analysis suggests that bacteria accumulate in fa- 
vorable regions if the impulse response function is nega- 
tive. As remarked in Ref. Q, this conclusion is paradox- 
ical in view of experimental measurements, which yield 
a bi-lobe response function Q. For comparable positive 
and negative lobes, chemotaxis ought to work best if the 
negative lobe is peaked around a time ro in the past, and 
fail if it is relegated much beyond in the past. We illus- 
trate this issue in Fig. 1, where we plot solutions of the 
non-linear model (Eq. ([4])) for a step of chemoattrac- 
tant concentration. We use a bi-lobe response function 
similar to the experimental one and derive the steady- 
state density of bacteria for three different values of ro. 
According to Fig. 1, accumulation in favorable regions 
occurs when ro is comparable to the time of the nega- 
tive peak in the response function (top curve in Fig. 1 
labeled by a disk symbol). For smaller values of ro, bac- 
teria feel the negative peak only rarely and accumulation 
occurs in unfavorable regions. This picture agrees with 
our analytical results in the linear model. 

Curiously, the experimental value ro generally quoted 
(~1 sec) falls between the two peaks of the response func- 
tion and, in our model, does not lead to favorable accu- 
mulation (intermediate curve in Fig. 1 labeled by an open 
square). This conclusion may be modified for a different 



shape of the response function, less similar to the exper- 
imental one — for example, one with a very deep negative 
lobe. Obviously, there are a number of constraints and 
performance requirements which we have not considered 
and which inform the shape of the single-bacterium fil- 
ter. For example, a rational for a response function that 
is spread out in time instead of narrowly peaked is the 
resulting robustness with respect to input noise, and a 
rational for a bi-lobe response function is the resulting 
'adaptive' mechanism of mean subtraction. 

In sum, we have introduced steady-state equations that 
govern bacterium density in chemotactic response to a 
chemoattractant profile. The solutions present a rich be- 
havior which depends in a subtle manner on the details of 
the model. We find that the bacterium density is a non- 
local function of the chemoattractant density (see Eqs. 
(HSHH]) and Fig. 1). This feature of the steady state is 
a direct consequence of the presence of memory in the 
dynamics and emerges in a proper treatment of correla- 
tions; earlier studies which ignore correlations find local 
solutions [1]. Our approach also predicts a regime in 
which bacteria accumulate favorably, even in the case of 
memory-less dynamics, in the tumbling state. Most ear- 
lier studies treat tumble as instantaneous. We treated 
tumble duration as a homogeneous Poisson process. In 
experiments, tumble duration seems to be influenced by 
the recent past in much the same way as run duration 
is, but with a bi-lobe response function that is more nar- 
rowly peaked and sign-inverted [3|. Roughly, we may 
say that tumbles tend to be shorter in favorable regions 
and longer in unfavorable regions. If so, chemotactic re- 
sponse may be weakened by this effect, with respect to 
the homogeneous tumble case. 
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